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Abstract 

We prove the bijectivity of the constraints of normalization and of the Fermi-Coulomb hole charge 
sum rule at each electron position for approximate wave functions. This bijectivity is surprising 
in light of the fact that normalization depends upon the probability of finding an electron at some 
position, whereas the Fermi-Coulomb hole sum rule depends on the probability of two electrons 
staying apart because of correlations due to the Pauli exclusion principle and Coulomb repulsion. 
We further demonstrate the bijectivity of these sum rules by example. 
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Sum rules play an important role in physics, and there are many ways in which they 
are employed. Within the realm of electronic structure theory, for example, accurate 
properties of a system may be determined by the use of an approximate potential whose 
parameters are adjusted so as to ensure the satisfaction of a sum rule. Metal surface 
properties such as the surface energy and work function are obtained by application of 
the Theophilou-Budd-Vannimenus sum rule jl| which relates the value of the electrostatic 
potential at the surface to the known bulk properties of the metal. The parameters in a 
model effective potential at a metal surface are then adjusted P] so as to satisfy this sum 
rule. Another manner in which sum rules have proved to be significant is in the context 
of Kohn-Sham density functional theory (KS-DFT) |3||, a local effective potential theory 
of electronic structure that is extensively employed in atomic, molecular, and condensed 
matter physics. In KS-DFT, all the many-body effects are incorporated in the 'exchange- 
correlation' energy functional of the ground state density. Since this functional is unknown, 
it must be approximated. A successful approach jj| to the construction of approximate 
'exchange-correlation' energy functionals, and of their derivatives which represent the local 
effective potential in the theory, is the requirement of satisfaction of various scaling laws 
together with those of sum rules on the Fermi and Coulomb hole charge distributions 
In the recently developed Quantal density functional theory (Q-DFT) |6|], the local effective 
potential is described instead in terms of the system wave function. Thus, one method for 
the construction of the local effective potential in Q-DFT is to employ an approximate wave 
function that is a functional of some functions Q. These latter functions are determined 
such that the wave function functional satisfies various constraints such as normalization, 
the Fermi-Coulomb or Coulomb hole sum rules, or reproduces a physical observable of 
interest such as the density, diamagnetic susceptibility, nuclear magnetic constant, etc. [3]. 

The satisfaction of a particular sum rule by an approximate potential, or an 'exchange- 
correlation' energy functional, or a wave function functional, however, does not necessarily 
imply the satisfaction of other sum rules. In this paper we describe a counter intuitive 
bijective relationship between the sum rules of normalization and that of the Fermi-Coulomb 
or Coulomb hole charge. The satisfaction of either one of the sum rules by an approximate 
wave function ensures the satisfaction of the other. This bijectivity is counter intuitive 
because the constraints of normalization and of the Fermi-Coulomb hole depend on 
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distinctly different quantum-mechanical probabilities. The bijectivity is also of importance 
from a practical numerical perspective. The proof and demonstration of the bijectivity of 
these sum rules constitutes the paper. 

To understand why this bijectivity is so counter to intuition, let us consider the physics 
underlying the two properties of an electronic system that these sum rules depend upon. 
For a system of N electrons, the constraint of normalization on an approximate wave 
function \I/(X) requires that 

y *(X)**(X)dX = 1, (1) 

where X = x 1; x N ; dX = dx 1; dx N ; x = r, s with r and s being the spatial and spin 
coordinates of an electron. (Atomic units e = Ti = m = 1 are assumed.) Equivalently, this 
sum rule may be written in terms of the electronic density p(r). The density p(r) is iV times 
the probability of an electron being at r : 

p(r) = NJ2 [ r^X^^^X^ 1 )^" 1 , (2) 
i J 

where dX 1 ^^ 1 = rfx 2 , dx N . The normalization sum rule then becomes 

J p{r)dT = N. (3) 

The density p(r) is a static or local charge distribution. By this is meant that its 
structure remains unchanged as a function of electron position r. Integration of this charge 
distribution — the normalization sum rule — then gives the number iV of electrons. Thus, 
normalization is a statement as to the number of electrons in the system. 

The definition of the Fermi-Coulomb hole charge distribution p xc (rr') derives from that 
of the pair-correlation density g(rr'). The pair-correlation density is the density at r' for an 
electron at r. The density at r' differs from that at r because of electron correlations due to 
the Pauli exclusion principle and Coulomb repulsion. Thus, the pair density is defined as 

g(rr>) = (^5^ - r)5(r, - r')|*>/p(r). (4) 
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Its total charge, for each electron position r, is therefore 

J g{rr')dv' = N - 1. (5) 

The pair-correlation density g(rr') is a dynamic or nonlocal charge distribution in that 
its structure changes as a function of electron position for nonuniform electron density 
systems. If there were no electron correlations, the density at r' would be p(r'). Hence, the 
pair-correlation density is the density at r' plus the reduction in density at r' due to the 
electron correlations. The reduction in density about an electron which occurs as a result 
of the Pauli exclusion principle and Coulomb repulsion is the Fermi-Coulomb hole charge 
p xc (rr'). Thus, the Fermi-Coulomb hole is defined as 

Pxc {vv')=g{vv')-p{v'). (6) 

The Fermi-Coulomb hole p xc (rr') about an electron is also a dynamic or nonlocal charge 
distribution. For nonuniform electron gas systems, its structure is different for each electron 
position. Since each electron digs a hole in the inhomogeneous sea of electrons equal in charge 
to that of a proton, it follows that the total charge of the Fermi-Coulomb hole surrounding 
an electron, for each electron position r, is 

J Pxc {vv')dv' = -1. (7) 

This is the Fermi-Coulomb hole sum rule. 



The definition of the Coulomb hole p c (rr'), which is the reduction in density at r' for an 
electron at r' because of Coulomb repulsion, in turn derives from that of the Fermi-Coulomb 
p xc (rr') and Fermi p x (rr') holes. The Fermi hole is the reduction in density at r' for an 
electron at r that occurs due to the Pauli exclusion principle. The Fermi hole is defined via 
the pair-correlation density g s (rr') derived through a normalized Slater determinant ${y?i} 
of single particle orbitals y?j(x): 

„ (rr H (^}lE^^-r)^-rO|${^}) 

9s[ } ~ p(r) W 

= p(v')+p x (vv'). (9) 

The orbitals <Pi(x) may be generated either through KS-DFT or Q-DFT in which case the 
density p(r) is the same as that of the interacting system, or they could be the Hartree-Fock 



theory orbitals for which the density is different. As the sum rule on g s (rr') is the same as 
in Eq. (5), and the Slater determinant is normalized, the total charge of the Fermi hole, for 
each electron position r, is also that of a proton: 



J Px {r,r')dv' = -1. 



(10) 



The Coulomb hole is then defined as the difference between the Fermi-Coulomb and Fermi 
holes: 

p c (rr') = p xc (rr')-p x ( r r f ). (11) 
. The total charge of the Coulomb hole, for each electron position r, is therefore zero: 



/ Pc {w')dv' 



0. 



(12) 



This is the Coulomb hole sum rule. 



Both the normalization and the Fermi-Coulomb or Coulomb hole constraints are charge 
conservation sum rules. However, their physical origin, and therefore the charge conserved 
in each case, is different. That these distinctly different charge conservation rules are 
intrinsically linked bijectively constitutes the theorem we prove. 



Theorem: The normalization and Fermi-Coulomb or Coulomb hole sum rules are bijec- 
tive. Satisfaction of the normalization sum rule by an approximate wave function implies the 
automatic satisfaction of the Fermi-Coulomb or Coulomb hole sum rules for each electron 
position. Conversely, the satisfaction of the Fermi-Coulomb or Coulomb hole sum rules for 
each electron position by an approximate wave function implies the normalization of that 
wave function: 

/ J^(X)*^(X)rfX= 1 \ / j p xc {rr')dr' = -1 \ 



or 

fp(r)dr = N 



or 



(13) 



^ / p c {vv')dv' = j 



Proof: (a) The proof of the arrow to the right in Eq. (13) is as follows. Let us assume 
an approximate wave function that is normalized. Then, integration of Eq.(6) over r' using 
the normalization constraint of Eq.(3) leads directly to the Fermi-Coulomb hole sum rule 
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of Eq.(7). 

(b) For the arrow to the left, consider an approximate wave function that satisfies the 
Fermi-Coulomb hole sum rule Eq.(7) for each electron position r. The sum rule Eq.(5) on 
the pair-correlation density g(rr') follows from its definition Eq.(4) which is independent 
of whether or not the wave function is normalized. Thus, since both the sum rules on the 
Fermi-Coulomb hole and the pair-correlation density are satisfied, then on integration of 
Eq.(6) over r', normalization of the wave function is ensured. 

(c) Consider an approximate wave function from which one constructs a Fermi-Coulomb hole 
for each electron position r. For a normalizd Slater determinant ${</?j}, next define a Fermi 
hole p x (rr') which then satisfies the Fermi hole sum rule of Eq.(lO). If the satisfaction of the 
Coulomb hole sum rule is now ensured, then this guarantees the satisfaction of the Fermi- 
Coulomb hole sum rule, which as shown in (b), ensures that the wave function is normalized. 

Recall that normalization depends upon the probability of finding an electron at 
some position. On the other hand, the Fermi-Coulomb and Coulomb hole sum rules 
depend on the reduction in probability of two electrons approaching each other. The 
fact that satisfaction of the integral condition of either one of these probabilities means 
the satisfaction on the integral condition of the other is not obvious, and therefore surprising. 

We next demonstrate the bijectivity of Eq. (13) by application to the ground state of the 
Helium atom. The nonrelativistic Hamiltonian of the atom is 

~ 1 o 1 o Z Z 1 . . 

H = -^ --V + — , 14 

2 2 n r 2 ri2 

where ri, r 2 are the coordinates of the two electrons, r 12 is the distance between them, 
and Z = 2 is the atomic number. The equivalence from left to right of Eq. (13) can 
be easily demonstrated by assuming an approximate wave function tp with parameters 
Ci(i = l,...,p) that is normalized in the standard manner at the energy minimized values 
of the parameters: = 0, where I[ip] = f ip*Hipdr/ f ip*ipdr. On the other hand, the 

equivalence from right to left is not as readily accomplished through such a wave function 
since the Fermi-Coulomb hole sum rule must be satisfied at each electron position. It is, 
however, possible to demonstrate the bijectivity by assuming the wave function to be a 
functional of a set of functions x- i> — i ) \x\ instead of simply a function. The functions 



6 



TABLE I: The satisfaction of the Coulomb hole sum rule Eq.(12) for different electron positions 



r(a.u.) 


/ p c (rr')dr' 


n 00566798 


-0 00039251 


0.13567807 


0.00032610 


0.57016010 


0.00034060 


0.72285115 


0.00013025 


0.89208965 


0.00001584 


1.07722084 


0.00007529 


1.49223766 


0.00029097 


1.96148536 


0.00034743 


3.91996382 


0.00032567 


5.15549169 


0.00057862 



X are determined so as to satisfy the normalization or Fermi-Coulomb hole sum rules as 
described in Ref.7. 

For the left to right equivalence, we choose the wave function functional to be of the form 

a 

i>[x} = ®(a,s)[l-f( X ;s,u)], (15) 

with$[a, s] = (a 3 /ir)e~ as , f(s,u) = e~ gu (l + qu)[l — x(l',s,u)(l + u/2)}, where a and 
q are variational parameters, s = r\ + r 2 ,w = ri 2 . The function x — X2 of [3], with 
the energy minimized values of the parameters being a = 1.6629, q = 0.17049. This 
wave function is normalized to unity, the function x being determined as a solution to 
a quadratic equation. We further assume, as in local effective potential energy theory, 
that the Fermi hole p x (rr') = — p(r')/2. The corresponding Coulomb holes p c (rr') are 
plotted in Figs. 1 and 2 for electron positions at r = 0,0.566,0.8, 1.0 (a.u.) together with 
the exact Coulomb hole 0. (The electron is on the z axis corresponding to 6 = 0. The 
cross section through the Coulomb hole plotted corresponds to 9' = with respect to the 
electron-nucleus direction. The graph for r' < corresponds to the structure for 9' = n and 
r' > 0.) The two Coulomb holes, though similar are inequivalent. Integration of both the ex- 
act and approximate Coulomb holes for each electron position leads to a total charge of zero. 

For the right to left equivalence of Eq. (13), we choose a different wave function functional 
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FIG. 1: Cross-section through the Coulomb holes for electron positions at (a)r = (a.u.), and (b) 
r = 0.566 (a.u.). The holes determined by the wave function functional of Eq.(15) and the 'exact' 
hole are plotted . 
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V4x] = -/(rir 2 )], (16) 

with /(nr 2 ) = e-^ 2 [l - X (R)(1 + r/2))] , r = n - r 2 , R = (n + r 2 )/2^ (3 = g[p(R)] 1/3 , q 
a variational parameter, and <&((j>i) the Hartree-Fock theory prefactor . The satisfaction 
of the Coulomb hole sum rule requires the solution of a nonlinear integral Fredho 
equation of the first kind for the determination of the function x(R). We have solved 
the linearized version of this integral equation for x(R). The satisfaction of the Coulomb 
hole sum rule for typical electron positions for q — 1 is given in Table I. (We do not 
plot the corresponding Coulomb holes as they are very similar to those of Figs. 1 and 
2.) The wave function functional of Eq. (16) thus determined satisfies the normalization 
constraint to the same degree of accuracy as that of the sum rule given in Table I. Hence, 
the bijectivity of the normalization and Coulomb hole sum rules is demonstrated by example. 
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FIG. 2: The figure caption is the same as in Fig.l except that the cross sections plotted are for 
electron positions at (c)r = 0.8 (a.u.), and (d) r = 1.0 (a.u.). . 



In conclusion, we have proved the bijectivity of the normalization and Fermi- Coulomb or 
Coulomb hole sum rules for approximate wave functions. The bijectivity is also significant 
from a numerical perspective because it is much easier to normalize a wave function than to 
ensure the satisfaction of the Fermi-Coulomb or Coulomb hole sum rules for each electron 
position. As shown by the examples, the determination of a wave function functional via 
normalization requires the solution of a quadratic equation, whereas that determined via 
satisfaction of the Coulomb hole sum rule requires the solution of an integral equation. On 
the other hand we note that the wave function functionals, as determined by satisfaction of 
the different sum rules, are different. Hence, the Fermi-Coulomb and Coulomb holes, and 
therefore how the electrons are correlated, will be different depending upon which sum rule 
is satisfied. It is unclear as to whether a better representation of the electron correlations 
is achieved by satisfaction of the normalization sum rule or that of the Fermi-Coulomb 
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hole. Finally, the bijectivity explains the results of our analysis ll| of the Colle-Salvetti 



wave function functional [12|. This wave function, which constitutes the basis for the most 
extensively used correlation energy functional in the literature, is of the same form as that 
of Eq. (16) except that x(R) = y/^P/iX + V^p), (3 = q[p HF (R)} 1 / 3 . In analyzing this 
wave function we had noted that it was neither normalized nor did it satisfy the Coulomb 
hole sum rule. These facts are consistent with the bijectivity theorem proved above. The 



lack of satisfaction of either one of the constraints ensures the lack of satisfaction of the other. 
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